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Abstract 

A type of super-hydrophobic surface consists of a solid plane boundary 
with an array of grooves which, due to the effect of surface tension, prevent 
a complete wetting of the wall. The effect is greatest when the grooves 
are aligned with the flow. The pressure difference between the liquid and 
the gas in the grooves causes a curvature of the liquid surface resisted 
by surface tension. The effects of this surface deformation are studied in 
this paper. The corrections to the effective slip length produced by the 
curvature are analyzed theoretically and a comparison with available data 
and related mathematical models is presented. 

1 Introduction 

Large pressure drops are necessary to cause liquid flow in micro- and nano- 
channels. The small values of naturally occurring slip lengths on hydrophobi- 
cally coated surfaces, typically of the order of some nanometers [1-4], are in 
most cases insufficient to significantly affect the required pressure gradients, as 
also confirmed by some theoretical analysis [5-7]. With the recent increasing 
interest in such systems [8,9], efforts to reduce these large pressures have been 
made. One promising avenue are the so-called super-hydrophobic surfaces [10- 
15], in which the wall stress is decreased by reducing the liquid-solid contact 
area [1,9]. The idea is to cover the surface of interest with structures such as 
posts, grooves, or others, over which the liquid surface remains suspended due 
to the effect of surface tension. The effect is equivalent to the introduction of an 
effective slip length of a magnitude comparable with the size of the geometrical 
features on the surface. This arrangement has been studied by several authors 
both experimentally [10-12,16] and theoretically [5,7, 17-22]. 
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Ou and co-workers [10,11] investigated the performance of a microchanncl 
with a rectangular cross section the floor of which contained a series of grooves 
aligned with the flow. This configuration resembles the situation studied the- 
oretically by Philip [19] who considered a variety of such arrangements. In a 
system of this type, the liquid surface becomes convex towards the gas occupying 
the groove due to the high pressure in the channel, with the radius of curva- 
ture determined by the balance of the pressure difference and surface tension. 
The importance of this effect, which was neglected by Philip, increases with 
increasing pressure and may be expected to become significant as the channel 
size decreases and the pressure levels correspondingly increase. The purpose of 
the present paper is to calculate the correction to the effective slip length due 
to this curvature. While the correction is small for the case considered by Ou 
and co-workers, one can easily envisage practical cases in which the effect would 
be important. Our approach is perturbative and rests on the smallness of the 
deformation of the free surface. 



2 Curvature effects on a periodic free shear pat- 
tern: formulation of the problem 

We consider a liquid flowing in a rectangular channel the floor of which contains 
equally spaced slots aligned with the flow. Surface tension prevents the liquid 
from filling the slots, but the pressure in the channel causes the liquid free 
surfaces to bow into them. If there is a sufficient number of slots, we can 
consider the flow as consisting of the periodic repetition of equal cells similar 
to the one shown in figure 1. In this figure, the top boundary, at a distance 
L* above the floor, is the roof of the channel, the dashed vertical lines are 
periodicity boundaries, and the lower boundary of the liquid region consists of 
a free-slip portion along a free surface S, occupying the range |ar*| < c*, and 
a no-slip portions on either side of it, for c* < |a;*| < H*/2, where H* is the 
dimensional cell width. Lengths will be rendered dimcnsionless with respect to 
H*/2ir, so that < c < tt and L = 2irL*/H*. It will also be convenient to 
introduce a slip fraction defined as 

(1) 

In fully developed conditions, only the axial velocity component w* is non- 
zero. Upon non-dimcnsionalizing it with respect to the axial pressure gradient 
\dp*/dz*\ and liquid viscosity (j,, w = (2n / H*) 2 /iw* / (— dp*/dz*), this quantity 
satisfies 

d xx w + d y yw = -1 , (2) 

where we have assumed the pressure field to be uniform over the cross section. 
The field w satisfies the no-slip condition on the roof of the cell and periodicity 
conditions on its lateral boundaries: 

w{x,L) = 0, d x w(±TT,y) = 0. (3) 
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On the floor of the cell, the no-slip condition applies away from the groove: 

w(x,0) = for c < \x\ < ir (4) 

while there is no tangential stress on the free surface S: 

n ■ Vw = on S (5) 

where n is the normal to the free surface. From a knowledge of w, one can 
calculate the flow rate in the channel 



Q 



[ w(x,y)dA (6) 

J A 



where A is the total cross-sectional area bounded by the solid walls and the free 
surface. 

One can define an effective slip length A by equating the actual flow rate Q 
with the flow rate Q e f / that would be found in a rectangular channel of height 
L and width 2n with a partial-slip condition applied at the bottom wall: 

Weff(x,0) = \d v W ef f{x,0). (7) 

The flow rate in this latter case is readily calculated: 

ttL 3 / A 



3 Perturbation problem 

With the small dimensions of practical interest, the curvature of the free surface 
is small and we represent it in the form 

y + e V (x) = (9) 

in which the parameter e, to be estimated presently, is taken to be small com- 
pared to 1. With this approximation it is easy to show that, correct to first 
order in e, the free-shear condition (5) is 

Tyz + T xz = (10) 

where t^ z = diW, i = x,y, denotes the only non-vanishing components of the 
viscous stress. The smallness of e suggests a perturbation approach to the 
solution of the problem and we write 

w = w (a) +ew (1) +o(e). (11) 

The field satisfies (2) while is harmonic. Both satisfy the conditions 
(3) and (4) on the solid and periodicity boundaries while, by a Taylor-series 
expansion in e, we find the boundary conditions in the groove as 

d y w {a) {x,0) = for \x\ < c (12) 
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and 



d y w {1) (x,o) = -a x ks)4» (0) (x,o) 



for 



\x\ < c 



(13) 



at orders e° and e respectively. When the expansion (11) is used to calculate 
the flow rate, we find 



+ Q 



(i) 



with 



and 



Q = Q (0) + e 
Q(°> = / da; / dyw^(x,y) 



o(e) 



Q{ 1] = [ dx [ dyw^(x,y) 

J-7T JO 

Q { 2 1} = J w(°\x,0)r){x)dx. 



(14) 

(15) 

(16) 
(17) 



This second term of the 0(e) contribution to the flow rate arises from the 
deformation of the free surface. 

In order to estimate the parameter e and the shape of the deformed free 
surface it is convenient to revert temporarily to dimensional quantities. From 
the pressure in the channel P* h and in the gas occupying the slot under the 
liquid, Pg as , we can calculate the radius of curvature R* of the interface from 
Laplace's formula as 



R* = 



a 



p* _ p* 

ch gas 



(18) 



where a is the surface-tension coefficient. The circle of radius R* passing 
through the points x* = ±c* , y* = has the equation 

2 



(x*) 2 + (y* - V(R*) 2 (c*) 2 ) = (R*f 



from which 



y 



(x*) 2 ). 



gas 



2a 



(19) 
(20) 

(21) 
(22) 



This relation can be written in the form (9) with 

. _ 1 
e 2R* 

so that 

IP 

47T 

Upon reverting to dimensionless quantities, the free surface shape r] in (9) 
is thus found to be given by 

r} = c 2 - x 2 . (23) 

The pressure in the channel falls in the direction of the flow, which will lead to 
an axial variation of R. This effect is usually sufficiently slow [10, 11] as not to 
significantly affect the assumption of parallel flow. 



4 



4 Zeroth order solution: Laminar flow over a 
flat patterned surface 

Following the approach proposed by Philip [19], we seek the solution in the form 

™ ( °' L) = -\v(v -L) + \lw^ l \x, y) (24) 

where the first term is the standard two-dimensional channel flow profile with 
no-slip top and bottom walls and w^°' L '(x, y) is the correction due to the free- 
slip portion of the lower boundary. Evidently u;(°< L ) satisfies Laplace's equation 
with conditions (3) and (4) on the solid and periodicity boundaries, while 

d y w(°' L) (x,0) = -1 for \x\ < c. (25) 

Since is periodic in x and even, it can be expanded in a Fourier cosine 

series in x, after which the requirement that it be harmonic determines the y 
dependence: 

(0,L) oc 

^(o,£)(z ; y) = ^o_ (1 _| ) + ^ a (o,L) [ x _ e -2n(L- w )] cos(nx) . (26) 

n— 1 

This form ensures the absence of slip on the channel roof y = L; the boundary 
conditions at the lower wall lead to the dual series equations 



(0,L) 



A0,L) 



J2 a n' L) (! - e- 2nL ) cos(nx) =0 c < x < tt (27) 



+ ^a( l ' i )n(l + e- 2nL )cos(nx) = 1 < x < c. (28) 

n=l 

This problem was solved by Philip [19], but his procedure was different from 
that adopted here and it is shown in appendix A that his solution is correctly 
recovered. The flow rate is readily calculated from (24) and (26): 

gm.^^^y (29) 

For a finite height L, it does not seem possible to solve the dual-series equations 
(27) and (28) exactly. We have calculated the solution numerically as explained 
in appendix C. The case of infinite depth can however be treated analytically 
(appendix A). In particular, one finds 

= -4 log (cos |) (30) 

and 

u; (0,oo) (x,0) = 2arccosh^^Jj^ \x\ < c. (31) 
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A quantity that is needed to compute the first order correction is d x w^°' L \x, 0) 
as it is evident from (13). This quantity can be easily expressed as 

oo 

d x wW {x, 0) = J2 na<£' L hm{nx) [e- 2nL - 1] (32) 

n=l 

and, in the limit L/H — > oo, the series can be summed exactly to obtain (see 
appendix A): 

sin — 

d x w^ oo) (x,0) = -V2^= 2 \x\<c. (33) 

^/cosx — cose 

5 First-order problem 

We seek the solution for the first-order correction to the velocity field in a form 
similar to (24): 

(1,L) oo 

w^ L \x,y) = ^-(1- j) + Y. ai n L) i 1 - e- 2n(L - v) ] e-^cos(nx). (34) 

According to (16), the contribution to the flow rate given by the velocity jW^' L ^ 
is 

Q^ = \a^L 2 . (35) 
The no-slip condition imposed on (34) is expressed by 

(0,L) oo 

^— + Y, ai n- L) ( 1 - e ~ 2nL ) cos ( nx )=° c<x<ir (36) 

2 n=l 

while the zero-shear condition (13) becomes 

(0,L) oo 

\r + ai n L)n i 1 + e ~ 2nL ) cos(nx) = d x {r)d xW (°' L \x, 0)) < x < c 

2L n=l 

(37) 

with d x w^°' L \x, 0) given by (32). Again, this dual-series system cannot be 
solved analytically except in the case of infinite depth for which we find (see 
appendix B) 

(i >oo) _ 4 [ c sm 2 (%)(<? -x 2 ) dx _ 2 f c (1 - cos(x))(c 2 - x 2 ) dx 



tt J cos(x) — cos(c) tt J cos(a;) — cos(c) 

(38) 

For general L, we have resort to a numerical solution. 
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6 Effective slip length 



We expand the slip length similarly to (14) as 



A (L) = X (0.L) + 



x (l,L) + A (1,L) +o(e) (3Q) 



in which A( 0,L ) is the contribution of the unperturbed flow (i.e., with a flat free- 
slip surface), A^ 1 '^ is the contribution of the perturbed velocity profile, and 
\^~' L ^ is the contribution of the unperturbed velocity arising from the deforma- 
tion of the flow passage. From (8), (14), (29) and (35) we find 

(0,L) 

A (0) = (4Q) 

X^ L ) =JL J C w W( x ,0)(c 2 -x 2 )dx. (42) 

Figure 2 is a graph of the zero-order slip length normalized by the width of the 
periodic cell: 

n(0,L) A (0,L*) 

= (43) 

2tt H* v ' 

as a function of the normalized channel height L/(2n) — L* / H* for various 
values of the slip fraction £. The straight lines are the corresponding analytical 
results for L/(2n) — > oo given by (see appendix A) 

: log cos — . (44) 



2tt tt ° V 2 , 

While convergence to this result is evident from the figure, it is also clear that the 
rate of convergence becomes slower and slower as the free-slip fraction increases. 
This feature is a consequence of the fact that, as follows from Eq. (7) and as is 
explicitly shown by (44), the slip length diverges to infinity when £ — ► 1. 

Figures 4 and 5 are a similar representation for the first-order corrections, 
A^' L ^ and X^'^, found from the numerical solution of (36) and (37), also for 
various values of the slip fraction £. Note that the term X^' L ^ (see figure 5) 
is a positive correction decreasing to zero when L/(2tt) — > oo. On the other 
hand, the negative correction A^ 1 is approaching a finite limit (straight lines 
of figure 4) in the limit of large channel height where, as shown in appendix B, 
le pa 1 

^ = -— F(£) 5* = ^- (45) 

2tt H* Kq) 2R* V ; 



e\i~'°°^ normalized to the pattern dimension can be represented as 
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with 

J a cos(stt^) - cos(tt£) 

A graph of this function is given in figure 3. Its asymptotic behaviors for £ near 
and 1 is given by 

f(o * y+£?+^?+o(e) e-o (47) 



and 



(48) 



7 Summary and conclusions 

Super-hydrophobic surfaces are necessary to significantly affect pressure gradi- 
ents and facilitate liquid flow in micron-scale channels. The super-hydrophobic 
effect is realized by patterning the channel walls with posts or grooves so that 
the liquid forms a free surface and remains partially suspended away from the 
wall due to the effect of surface tension. While, on small scales, surface tension 
is a powerful force, it cannot prevent the free surface from bulging out becoming 
convex toward the gas space. This circumstance has two effects. In the first 
place, the flow passage is enlarged, which increases the slip length but, on the 
other hand, the velocity field is modified, which decreases it. 

In this paper we have considered one special situation of this type - a flat wall 
with grooves parallel to the flow direction - thus providing a natural generaliza- 
tion of the analysis carried out by Philip [19] who assumed a flat the liquid-gas 
interface. We have quantified the two effects mentioned before finding that the 
magnitude of the first one relative to the second one decreases proportionally 
to the ratio of the pattern width to the channel height. For deep channels, the 
more significant effect is the second one, given in (45) and shown in figure 3, 
which we found to give a negative contribution to the slip length. Our results 
are described in detail in section 6. 

The magnitude of the effect that we have studied is quantified by the dimen- 
sionless ratio 

IT* D* D* 

e = — (49) 

4tt a 

The only existing data with which our theory can be compared are those re- 
ported by Ou and coworkers [10]. We can estimate from their figure 9 that the 
air-water interface protrudes 2 to 4 fim. below the channel floor. The spatial 
period of the grooves on their lower wall is H* = 60 /um. Our parameter e is 
therefore of the order of (2-4)/60 = 0.03-0.06, which produces only a very small 
correction to the unperturbed solution. However, for periodicity patterns H * ~ 
10 /im, with P* h - P* as ~ 100 kPa, a ~ 0.1 N/m, we have e ~ 0.1, an estimate 
that would increase further with increasing area fraction £ as shown in figure 3. 
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8 Appendix A 



In this appendix we solve the zeroth order problem for a channel of infinite 
depth. In this limit, the appropriate situation to consider is that of a linear 
shear flow over a periodic array of free-slip longitudinal strips. We write the 
velocity field as the sum of a linear shear flow plus a perturbation (W ' 00 '): 



w^ co \x,y)=y + w^ 00 \x,y) 
with the correction expanded as 



(0,00) 00 

t(°'°°\x,y) = %— + Y / a i a'° c) cos(nx)e- ny . 
2 „=i 



w 



The boundary conditions (4) and (5) give rise to the dual series problem 

(0,00) 00 

+ V 4 0lOO >cos(na;) - c < x < tt 
2 „=i 



^2 a^ ,00 ^ncos(nx) = 1 < x < c. 



71=1 



For the problem 



00 

^4 0,OO) + al 0,Oo) cos(nx) =0 c < x < n 



n=l 



na( l °' oo) cos(nx) = f(x) < x < c 



Sneddon (p. 161) [24] gives the solution 



(0,00) _ 2 



(50) 
(51) 

(52) 
(53) 

(54) 
(55) 

(56) 



,(0,00) 



h ( °'°°\t)[P n {cos(t)) + P n ^(cos(t))}dt 



.2V2J0 

with P n Legendre polynomials. The function /4°' 00 ' ) (£) is 



n= 1,2,.. 



sin(^) dx 



n dt J ^/cos(a;) — cos(i) Jo 



f(u)du 



(57) 
(58) 
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or, m our case, 

ft 



(0,00)^ 2 d f xsin(|) dx 



h?'°°>(t) = -^- - K2; =. (59) 
^ dt Jo y'cos(x) — cos(i) 

The computation of ao°' °' ) and of u^ 0,00 ) hinges on the knowledge of the 
function h^'°°\t) which is the derivative of 

rn~,\, s 2 /"' a; sin(§) da; , % 

7 (0 '° o) (*) = - / , ; ^ (60) 
n Jo \/coa(x) — cos(t) 

This integral can be evaluated with some manipulations and the use of formula 
3.842 of Gradshteyn & Ryzhik [23] with the result: 



from which 



and 



^°' oo) (t) = x/2tan^ (62) 



,(0,oo) / 1 



5— ^^l- (63 » 



We also notice that from Sneddon (p. 161) [24] 



) + Val°^cos(nx)=cos(J) f (64 ) 
2 ° ^ V 2;; x ^cos(x) - cos(t) 

that immediately leads to the velocity at y = 

t/)(°-~>(M) = V^cq S (^) f-p tan ^^ dt \x\<c. (65) 

J x y/ cos(x) — cos(i) 

This integral can be done exactly leading to (31). The quantity d x w^ '°°\x, 0) 
is also readily evaluated with the result given in (33). 

9 Appendix B 

In this appendix we solve the first order correction to the linear shear flow 
problem considered in appendix A. We write the velocity corrections as 

(l.oc) CO 

(x, y) = + V a^cos{nx)e- ny (66) 

2 n=l 

with the boundary conditions (4) for c < x < tt and (13) for < x < c. 
Formulae (56) and (57) of appendix A again apply with the function h\ now 
given by 

(li00) 2 d f< S m(%Mx)d x w< Q >°°\x,0) dx 

h \ v) = zt+ / ? s m ( 67 ) 

tt dt J a y/cos(x) - cos(f) 



-a 



10 



from which wc have 



a o / — 1^^^=^^^= 

n Jo ^/cos(x) — cos(c 



(68) 



or, upon using (33) for d x w^ 0, °°^ (x, 0), 

sin 2 (f)(c 2 -a; 2 ) dx 2 f c (1 - cos(x))(c 2 - x 2 ) dx 



o (i,oo) = _4 r sin z (|)(^-^) _ _2 /" 
' :i 7T 7 cos(a;) - cos(c) tt J 



ir J cos(x) — cos(c) 

(69) 

If we introduce the slip length as Ai 1 ' 00 -* = ^ — , when we express eA^ 1,00 ' 
normalized to the pattern dimension we obtain: 



A< 1 ' oo) 1 T (1-cos(x))(c 2 -.t 2 ) da; 



- -— f 



27r 47T 2 i? 7 cos(:r) — cos(c) 

If we introduce s = x/c, 5=-^ and c = tt£ the previous expression becomes 

A (l,oo) g 



(70) 



e- 



7T F (0 (71) 



2tt 2tt 
with 

F(n = e C (i-«>^0)(i-^ 2 ) r721 

A partial evaluation gives 

F ® = ^ tan fa) (j^ S ^ da ~ <Si( 2 <)) - ( 73 ) 

in which 

_ , . v-^ sin(na) 1 f a 1 , 

s,(a) = s "V - 5 y„ io8 5(1^)) % - (74 » 

Integration by parts leads to 

/ 1 / 1 \ \ f a ' 2 

S\(a) = a ^log - — log sin ^-aj J + 2 J xcotxdx. (75) 

The last integral cannot be evaluated in closed form. 

The asymptotic results mentioned at the end of section 6 are readily derived 
from these expressions. 



10 Appendix C 

In this section we briefly sketch the procedure used to solve the dual series equa- 
tions numerically. Both the unperturbed problem (27), (28) and the perturbed 
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problem (36), (37) have the general structure 

A 00 

— - + ^ ^nA n cos(nx) = c < x <tt (76) 
2 n =i 

OO 

J~] n^ n A n cos(nx) = f(x) 0<x<c (77) 

n=l 

with r„ and generic functions of n. The two expressions are multiplied by 
cos(n;r) and integrated in their respective domains of validity. Upon using the 
identity 

/ cos(mx)cos(nx)dx — -ir5 nm — / cos(ma;)cos(nx)rfa; (78) 
Jo 2 J c 

the result may be written as 

f V A 00 ^ / rc \ 

/ — cos(mx)dx + \ J r„A„ I — / cos(ma;)cos(na;)dx + — <5„ r „ I = 
Jc 2 n=1 \ Jo 2 y 

y rVS> n A n I cos(nx)cos(mx)dx — / cos(mx) f (x)dx 

Upon adding these two relations, the result may be written in the form of a 
linear system: 

00 

Y, M n ^ m A n = B m (81) 

n=0 

where 

1 r 

Mo, m = - / cos(ma;)dx (82) 

2 J c 

r I" 

M„, m = (n*„ - r„) / cos(mx)cos(nx)dx + -S nm T n (n > 1) (83) 

Jo 2 

B rn = / cos(mx)f(x)dx (84) 
Jo 

The linear system (81) has been truncated and reduced to a N x TV matrix and 
is then found to converge upon truncation refinement. 



(79) 
(80) 



12 



References 



[I] E. Lauga, M. Brenner and H. Stone, "Microfluidics: The no slip boundary 
condition" in Handbook of experimental fluid mechanics (Springer, 2006). 

[2] J. T. Cheng and N. Giordano, "Fluid flows through nanometer scale- 
channels", Phys. Rev. E 65, 0312061 (2002). 

[3] O. Vinogradova, G.E. Yabukov: "Dynamic Effects on Force Measurements. 
2. Lubrication and the Atomic Force Microscope" Langmuir 19, 1227-1234 
(2005). 

[4] Y. Zhu and S. Granick," Limits of the hydrodynamic no-slip boundary con- 
dition", Phys. Rev. Lett. 88, 106102 (2002). 

[5] J.-L. Barrat and L. Boquet ," Large slip effect at a non- wetting fluid-solid 
interface " Phys. Rev. Lett. 82, 4671-4674 (1999). 

[6] C. Cottin-Bizonne, C. Barentine, E. Charlaix, E. Boquet and J.-L. Barrat, 
"Dynamics of simple liquids at heterogeneous surfaces: Molecular dynamics 
simulations and hydrodynamic description" , Eur. Phys. Jour. E 15, 427-438 
(2004). 

[7] R. Bcnzi, L. Biferale, M. Sbragaglia, S. Succi and F. Toschi,"On the 
roughness-hydrophobicity coupling in micro- and nano-channel flows" , Phys. 
Rev. Lett, submitted (2006). 

[8] C.-M. Ho and Y.-C. Tai, "Micro-clcctro-mcchanical systems (MEMS) and 
fluid flows", Annu. Rev. Fluid. Mech. 30, 579-612 (1998). 

[9] P. Tabcling, Introduction a la microfluidique (Bclin, Paris, 2003). 

[10] J. Ou, B. Perot and J. Rothstein," Laminar Drag reduction in microchannels 
using ultra-hydrophobic surfaces", Phys. Fluids 16, 4635 (2004). 

[II] J. Ou and J. Rothstein," Direct velocity measurements of the flow past 
drag-reducing ultra-hydrophobic surfaces" Phys. Fluids 17, 103606 (2005). 

[12] P. Joseph, C. Cottin-Bizonne, J.-M. Bcnoit, C. Ybert, C. Journct, P. Tabel- 
ing and L. Bocquet," Slippage of water past superhydrophobic carbon nan- 
otubes forests in microchannels" Phys. Rev. Lett. 97, 156104 (2006). 

[13] J. Bico, C. Marzolin and D. Quere," Pearl drops" Europhys. Lett. 47, 220- 
226 (1999). 

[14] D. Oner and T.J. McCarthy," Ultra-hydrophobic surfaces. Effects of topog- 
raphy length scale on wettability" Langmuir 16, 7777-7782 (2000). 

[15] K. Watanabe, Y. Udagawa and H. Udagawa, "Drag reduction of Newtonian 
fluid in a circular pipe with a highly repellent wall" Jour. Fluid. Mech. 381, 
225-238 (1999). 



13 



[16] C.-H. Choi and C.-J. Kim, "Large Slip of Aqueous Liquid Flow over a Na- 
noengineered Supcrhydrophobic Surface" Phys. Rev. Lett. 96, 066001 (2006). 

[17] R. Benzi, L. Biferale, M. Sbragaglia, S. Succi and F. Toschi,"Mesoscopic 
modelling of heterogeneous boundary conditions for microchannel flows", 
Jour. Fluid. Meek. 548, 257-280 (2006). 

[18] E. Lauga and H. Stone, "Effective slip length in pressure driven Stokes 
flow", Jour. Fluid Mech. 489, 55-77 (2003). 

[19] J. Philip, "Flow satisfying mixed no-slip and no-shear conditions", Z. 
Angew. Math. Phys. 23, 353-370 (1972). 

[20] M. Sbragaglia and A. Prosperetti," Effective velocity boundary condi- 
tion at a mixed slip surface " Jour. Fluid Mech., submitted (2006). See 
physics/0607003. 

[21] N. V. Priczjcv, A. A. Darhuber and S. M. Troian "Slip behavior in liquid 
films on surfaces of patterned wettability: Comparison between continuum 
and molecular dynamics simulations" Phys. Rev. E71, 041608 (2005). 

[22] J. Davies, D. Maynes, B. W. Webb and B. Woolford," Laminar Flow in a 
microchannel with superhydrophobic walls exhibiting transverse ribs" Physics 
of fluids 18, 087110 (2006). 

[23] I. S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, 
6th edn. San Diego: Academic Press (2000). 

[24] I.N. Sneddon, Mixed boundary value problems in potential theory (Amster- 
dam: North-Holland, 1966). 



14 




Figure 1: Schematic geometry of the problem studied in this paper. The floor 
of a channel with height L* contains a series of regularly spaced grooves parallel 
to the flow direction. When the number of grooves is large, the situation can 
be approximated by the periodic repetition of a fundamental cell of width H* 
such as the one shown here. The width of the the groove, where the shear stress 
is essentially zero, is 2c*. For small deformation, the maximum penetration of 
the free surface S into the groove is S* = (c*) 2 /2R* where R* is the radius of 
curvature of the free surface. 
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Figure 2: The zero-order effective slip length A^ '^ normalized by the horizontal 
period 2tt as a function of the channel height L also normalized by 27r for various 
values of the slip fraction £ = c/ir; the nomenclature is defined in figure 1. The 
symbols are results obtained by solving numerically (27) and (28). The straight 
lines are the limit L/(2n) — > oo given by the exact solution (44). 
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Figure 4: Results for the quantity A^ L '/(27r) as a function of the dimensionless 
ratio L/(2it) for various values of £. The symbols are results obtained by a 
numerical solution of (36) and (37). The limit L/(2n) 3> 1 corresponds to the 
analytical solution computed in appendix B. 
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Figure 5: Results for the quantity X^'^ /(27r) as a function of the dimensionless 
ratio L/ (2tt) for various values of £: £ = 0.5 (top) £ = 0.45 (middle) and £ = 0.4 
(bottom). Results are obtained by a numerical solution of (36) and (37). 
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